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Abstract 

Considering corrections from two-loop Feynman diagrams which involve gluino at large tan f3, 
we analyze the effects of possible CP phases on the rare B decays: — > and B Kl~^l~ in 
the CP violating minimal supersymmetric extension of the standard model. It is shown that the 
results of exact two loop calculations obviously differ from that including one-loop contributions 
plus threshold radiative corrections. The numerical analysis indicates that the possibly large CP 
phases strongly affect the theoretical estimation of the branching ratios, and this results coincide 
with the conclusion of some other works appearing in recent literature. 
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I. INTRODUCTION 



The rare B decays serve as a good probe for the new physics beyond the standard 
model (SM) since they do not seriously suffer from the uncertainties caused by the long 
distance effects. The forthcoming experiments of the B factories will make more precise 
measurements on the rare processes and it is believed that those measurements should set 
even stricter constraints on the parameter space of new physics. Among those plausible new 
physics scenarios beyond the SM, the simplest and most favorable extension is the Minimal 
Supersymmetric Model (MSSM) jl|. In the MSSM, there are five physical scalars (Higgs) 
compared to the SM where there is only one. The contribution of those neutral Higgs bosons 
to the rare B processes have been extensively discussed in literature. The main conclusion 
of the analysis is that the branching ratios of the processes such as B^ — > B^ — > Kl^l~ 
etc. are enhanced for larger values of tan (3, which is the ratio of the two vacuum expectation 
values (VEVs) of the neutral Higgs fields, even within the minimal flavor violating (MFV)^ 
supersymmetry scenario. 

In addition to the SM CP phase which exists in the CKM matrix, there are three more 
possible sources for the CP violation phases in the MSSM Lagrangian. The first one is the ix 
parameter in the superpotential which is complex and the second source is the corresponding 
soft breaking parameters, namely, the complex masses of the SU (3) x SU (2) x U{1) gauginos 
in the soft breaking terms induce three CP phases. As the third one, there are several CP 
phases emerging from the scalar soft mass matrices u d l r trilinear coupling 

matrices A^^^. Generally, only the off-diagonal elements of the soft mass matrices can be 
complex due to the hermiticity of these matrices. By contrast, the trilinear coupling matrices 
Ay ^ ^ can have complex diagonal elements j2]- Actually, the CP phases must be constrained 
by the present experimental results. The most rigorous constraints on those CP phases come 
from the experimental bounds electron and neutron electric dipole moments (EDMs), 

which are < 4.3 x 10~^^e ■ cm [3] and < 6.5 x 10~^^e ■ cmQl, respectively. The bound of 
the EDM of H^^^: d^^^^ < 9. x 10~^^e ■ cm[5j is also measured with high accuracy. In order 
to make the theoretical prediction consistent with the experimental data, three approaches 
are adopted in literature. One possibility is to make the CP phases sufficiently small, i.e. 



^ That is, the Cabibbo-Kobayashi-Maskawa (CKM) matrix is assumed to be the only source of flavor mixing. 
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n 

smaller than 10 ^ [6j. Alternatively one can also assume a mass suppression by making 
the supersymmetry spectra heavy enough, i.e. in a range of several Tev[3], or invoke a 
cancellation among different contributions to the fermion EDMs 
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12| | analyzes the mixing among the 



In the second scenario, a series of works[ 
neutral Higgs bosons in the CP violating MSSM. Considering the constraints from the 
experimental upper bounds of the electron and neutron electric dipole moments (EDMs), 
the soft trilinear coupling for the third generation scalar quarks can have large CP phases. 
Due to the large Yukawa couplings for the third generation quarks, the radiative corrections 
can lead to a large mixing among the would-be CP-even and CP-odd neutral Higgs bosons. 
This mixing causes drastic changes for the couplings between the neutral Higgs and quarks, 
the neutral Higgs and gauge bosons, and the self couplings among the Higgs fields, thus the 
consequence is that the lower bound of the mass of the lightest neutral Higgs is relaxed to 
60 GeV. 



Presently, the calculations of the effective 
presented by the authors of Ref . 3: 
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lamiltonian for the transition b sl~^l are 
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18| within the MSSM and THDM models. 



However, those works are all focusing on the CP-conserving processes. As pointed out in 
R,ef . [ifljl . the CP violation phases which induce mixing among the neutral Higgs bosons 
affect the effective Hamiltonian of 6 ^ sl^l~ significantly, especially for larger tan (3 values. 
However, their analysis only included the leading terms in tan^ (3 which originate from 
the counter diagrams. 

In this work, we analyze the rare S-decays: B X^'-f, l^l^ and B Kl^l~ 

{I = /i, r) in the CP violating MSSM. In order to reduce the number of free parameters, 
we assume no additional sources of flavor violation other than the CKM matrix elements. 
In our calculation, we consider all possible contributions from the one-loop box, 7—, Z— 
and Higgs-penguin diagrams, as well as the threshold radiative corrections to the effective 
Hamiltonian for b sl^l^ at first. The threshold radiative corrections indicate a sum of the 
SUSY QCD and SUSY electroweak corrections to the effective vertex H*D'^Q at larger tan/? 



values j20|, |21|, 122, |23[. Although threshold radiative corrections approximate the exact two- 



loop results adequately when the supersymmetry energy scale is sufficiently high, the authors 



of Ref. 



n 



2J, pointed out that the difference between threshold radiative approximation 
and exact two loop calculation is obvious in some regions of the parameter space of the 
MSSM. As a comparison and complement to the results of threshold radiative corrections 
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which were aheady derived and evaluated in hterature, we present calculations of the exact 
two-loop corrections from gluino at large tan (3 here. Those two-loop results, which include 
the gluino corrections to the penguin vertices shZ , sbH, and the gluino corrections to the 
four fermion interaction sbll, are explicitly presented in this work. The gluino corrections 
to the penguin vertex sb'j can be found in our previous analysis j^]. 

In the second section, we will present the modified couplings involving the neutral Higgs 
bosons in the CP violating MSSM. The effective Lagrangian for rare B decay modes is given 
in the third section. We also show the newly derived theoretical formulations of those decay 
branching ratios etc. in this section. Section IV is devoted to the numerical analysis and 
discussion, our conclusion is made in the last section. Some tedious formulae are collected 
in the appendices. 



II. MIXING AMONG NEUTRAL HIGGS BOSONS AND RELEVANT COU- 
PLINGS 

In this section, we present the modified Yukawa couplings of the neutral Higgs to the 
squarks and sleptons where the effects induced by the CP violation phases are taken into 
account. The most general form of the superpotential which has the gauge invariance and 
retains all the conservation laws of the SM is written as 

where are the Higgs superfields; and are quark and lepton superfields in 

doublets of the weak SU(2), where 1=1, 2, 3 are the indices of generations; the rest superfields 
?7^, t)^ and are the scalar quark superfields of u- and d-types and charged leptons in 
singlets of the weak SU(2) respectively. Indices i, j are contracted for SU(2), and h\ /i"''^ are 
the Yukawa couplings. To explicitly break supersymmetry, the soft-supersymmetry breaking 
terms are introduced as 

C^^f^ = -m^ ^ Hi* Hi - ^ H\*H1 - ^ LI*LI - ^ R'*R' - ^ Ql*Ql - ^ U'*U' 

-m'^^^D^*D^ + (miAsAi + m2X\X\ + m^X^X^ + h.c.) + [j^Be^jHlH] + eijA^HmPi^ 
+e,jA''HlQ]D' + ei,A^H^Q]u' + h.c] , (2) 

where m? , m? , m? , , m? , and are the squared masses of the superpar- 
tides, ms, m2, mi denote the masses of (a = 1, 2, ■■■ 8), A^ [i = 1, 2, 3) and Xb, 



which are the SU{3) x SU{2) x f/(l) gauginos. 5 is a free parameter in unit of mass. With 
the soft-supersymmetry breaking terms in Eq.Q, we can study the phenomenology in the 
framework of the MSSM. 

By the effective potential which accounts for the two-loop Yukawa and QCD corrections 
via the renormalization group equation (RGE), the squared mass matrix of the neutral Higgs 
bosons is written as: 



m 



( ( 



2 2 

a ^ 



+ \s„c 



hRe{A6)c^ + Re(A7)s^ 



\ Im(A5)s^ + Im(Ag)c^ 

with the squared mass being 



2 2 

m„ = m , 



^[i + A4 s,, c 



13^ U 



+ Re(Ae)c^^ + Re(A7)sJ^ 



Im(A5)c^ + Im{A7)s^ 

4m^s^ /I 
^2 



Im(A5)s^ 
+ Im(A6)c^ 



Iiii(A5)s 



|3 



+ Im(A6)c^ 



(3) 



(4) 



where the parameters represent the masses of the physical charged Higgs-bosons. The 
concrete expressions of the parameters Aj {i = 1, 2, ■ ■ ■ , 7) can be found in Ref. jl^. Since 
the squared mass matrix of neutral Higgs m^^ is symmetric, we can find an orthogonal 
matrix to diagonalize it: 



Z m Z^ = diaqim ,m ,m ) . (5) 
Correspondingly, the modified interactions related to our calculation are listed below: 
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em 
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em. 



C 



e 



■(C) .i^i 



2i 



0~ ~* 



s^Z^' + c,Z^' + iZ'' 

13 H 13 H H 



li 



s.Z''' + c.Z''' - iZ 

13 H 13 H H 



W;{H'Sd'H-)-{id^H'',)H- 
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em 



"0 w=l= w"F 



H P H 



(6) 



The couplings k^'^, C'' are presented in appendix O Since the SUSY-QCD modifies 



the Yukawa couphng tbH^ remarkably 



22j | for larger tan/5 values, we should consider the 



SUSY-QCD effects on the effective lagrangian which only includes one loop contributions. 
Re-summing the dominant supersymmetric corrections for larger tan (3 to all orders, as well 
as the leading and next-to-leading logarithms of the standard QCD corrections, one can 
write the interaction as 2^, 0| 

em,(g) 



C 



Hqq 



2m s s„c„ 

vj (3 (3 



tb'^FC 



3k 



Hi. 



e 



+ AZf + - c^A)zl^]b^~bHl + h.c. 



^ m 

w 



* ^ -buj+t + ^1 "^^^ tan I3buj-t 



m„ tan /? 



m„ 



v: 



buo^t 



m„ 



buj+Ui H tan pbu^Ui 



V* 

CKM ib 



buj+Ui — 



tan P 



m 

w 



H- 



m„ 



m„ 



* ^CKM ti 



+ 



in.n(3diU0_i\H- + y C* V* 



m,. 



m^tan/3 



where Q is the characteristic energy scale of the process. The tan /3-enhancing radiative 
corrections are 

1 

X 



X 



l + tan/3A' 

em^{Q) tan f3A^ 



FC 



-X, 



v^m„s„s^[l + tan/5A 
e = xJl + A'-/i, cot/3A" 



C = xJi + ^ 



tub ^-cfe 

^CKM ^CKM 



^CKM ^CKM 



1 + tan /3A 
1 + tan /JA"" 



cud 
'CKM 



A A S , . EW 

A = A + /l,A 



A^ 



2c^ 



m*a*I(m. ,m. .ImA) 



647r2 



(8) 



with the vertex function 
/(a, 6, c) 



a2-62)(62-c2)(a2-c2) 



Vln(^)+6Vln(^)+cVln(^)]. (9) 



While deriving Eq. ((Tj) and Eq. (jHJ, we consider the fact that ^ (m^, m^); S> 
(m^, m^) and |V^J (|V^J, |V^^|, |V^J, 1Kb I)- Computations without these approximations 
were given m Ref.Q. The running quark masses are evaluated by 



(10) 



where we have assume that there are no other colored particles with masses between Q and 

aAQ2)\df 



m^. The evolution factor reads 



f/,(Q2,Qi 



d 



12 



J 



^ 33 - 2/ ' 

8982 - 504/ + 40/2 



(11) 



3(33 -/)2 

where, / is the number of active quark flavors. With those preparations given above, we 
can discuss the rare B decays in the CP violating MSSM. 



III. THE EFFECTIVE HAMILTON AND DECAY WIDTH FOR RARE B DE- 
CAYS 

A. The effective Hamiltonian 

The processes which we are interested in, are —>■ l^l^ and B Kl^l^, both of 
them originate from the transition b ^ s. Integrating out the heavy degrees of freedom in 
the full theory, an effective Hamiltonian is obtained 26[: 

4C I- 10 10 . 
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with Coim ) = ff*^ , v.. represent the physical CKM entries, and f** is the 

tan /5— enhanced radiative corrections to the effective CKM entries. The operators in the 
effective Hamihonian are: 



01- - 

02- - 

03- - 

0-J-- 
Os-- 

O'- 



O' 
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o' 



O' 



^ oi~1 —^a 

em 



J2 [Qal^^-Qis 

q=u,d,s,c,b 
q=u,d,s,c,b 
q=u,d,s,c,b 
q=u,d,s,c,b 



(4vr) 
em 



(4vr) 



2 "'a-' 0/3 



+ "/3 , 



: ^(s„7M^+&a)(/"7'^75/ 

^(s,^+6,)(r75/) , 

g2 
g2 

^(.-„^_6„)(r75/) . 



(13) 



In our calculations, we adopt the Feynman rules in the 't Hooft-Feynman gauge with 
^ = 1. Within the framework of the MFV CP violating MSSM, the one loop Feynman 
diagrams that contribute to the effective Hamiltonian (Eqll2|l can be found in the literature. 
For example, the one- loop 7—, Z— penguin diagrams are presented in Ref. and the 
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'z, G°, H° 

(a) 



; Z, G°, H° 
(b) 




FIG. 1: The corrections of two-loop diagrams involving gluino (later we abbreviate them as "the 
two loop gluino corrections" ) to the penguin vertices sbZ, sbH at large tan /?. 

Higgs penguin and box diagrams are given in Ref.j3|. At the electro- weak scale, the one 
loop Wilson coefficients in Eq^Jare divided into several pieces: 

= CUf^^) + Cf{f,J + {i = 9, 10) 

ClM = + Cf (/"w) + a'""(/iw), = 9, 10) 

= cf^M + cr'^'M + c^^(/iw) + c^r^^™(/iw), = s, P) 
ciM = c'rM + cr"*(/^w) + c':>"^M, {^ = s,p). (m) 
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Here, we collect the Wilson coefficients corresponding to the one loop 7—, Z— , H^— pen- 
guin and box contributions in the appendix El In order to include the threshold radiative 
corrections of gluinos, we have replaced the tree-level vertices by the corresponding inter- 
actions which are presented in Eq. (|7|). Similarly, we can obtain the contributions which 
originate from a resummation of high order threshold effects on the Wilson coefficients of 
the operators involving down-type quarks: 



C^-^u ) = -z^^x* Z'^U.^'" - s.Z'" + iZ^'^\ , (15) 

P VA^W/ 2C X PC H \ l3 H 13 H H J ^ ^ ' 



with = cos/3, = sin/3, and = 

w 

As we mentioned in the introduction, there are obvious differences between the results 
from exact two-loop calculations and that from threshold radiative approximation which 
is derived in terms of the heavy mass expanding method. In the large tan/3 scenario, 
corrections from the two-loop Feynman diagrams including gluino to the penguin vertices 
sbZ, shH are drawn in FigO Correspondingly, the corrections to the Wilson coefficients 
from the two-loop penguin diagrams are written as 

w w 

w w 

w w p — i p Ml-' 



I \ a "fnj^i^a OL m.m ^tl JL / x 1 



(16) 



8 ... 10 



with „ = 11 'P^^^ ^ 'P o = 'P^^^ ■ We list the expression of those nonzero form factors 

i=l " i=l h" 

-pC*) -pC*) in the appendix^ in order to shorten the length of text. 

Z.G fjP 

We should also include the two loop gluino correction to the four fermion interactions 
sbll for completeness. The Feynman diagrams are plotted in FigEl Correspondingly, the 
corrections to the Wilson coefficients from those two loop box diagrams are formulated as 

W 
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FIG. 2: The two loop gluino corrections to the four fermion interactions sbll at large tan/3. 



w 
w 
w 
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S C' 



with 



2a 

w 

-^t B' 



(17) 



B' ^ = Bf] + Bf] 



B^ 



B' 



^-^ S,P 

i=l 

6 

/-^ S.P 



(18) 



i=l 



The concrete expressions of those nonzero form factors SjJ, Bfl can be found in the 

appendix O In Eq. UHland Eq. [T71 we adopt the MS scheme to remove the UV-divergences 
which are caused by the divergent sub-diagrams 0]. 

When the effective Lagrangian is apphed to the hadronic processes whose characteristic 
energy scale is about m^, we should evolve those Wilson coefficients from the weak scale 
down to the hadronic scale. The running depends on the anomalous dimension matrix of 
the concerned operators [3(|. The Wilson coefficients obtained at the weak scale are regarded 



as the initial conditions for the renormalization group equations (RGEs). Up to the leading 
order (LO), the Wilson coefficients at hadronic scale are given as 



lu^s ■ u p 

aii,y,y] 



16 o 2021 

Crim,) = r/-i[C7(m^) + ^(Csim^) - -^^2(171^)) {v^^ - l) 

389^ , ^ / 10 \ 107^ , X / 28 -.^ 

+ ^C2{m^)[V^^ - 1) + 7;7^C'2("^w)(^'^ - 1 



540 



702 



Cg{mJ = C9(m^) + 
Cio("^J = Cio(m^) . 



An 



4 . lis 8 , 



29 , 
23 



(19) 



with T] = The most stringent constraint on the supersymmetry parameter space 

comes from the rare decay of 5— meson: B X j Isjl. The theoretical prediction on the 

nn 

branching ratio of the inclusive process B — > X^^y is given as |33L I35l| 

TiB ^ X 7) 



BRiB ^ X 7) 



TiB ^ X^eu^ 



-BR{B ^ X ez/^ 



\V*V P 



b 



BR{B X,ez7, 



(20) 
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where a^^ is the QED fine structure constant and the phase factor is 



/"ITL\'2 



m\6 /m 



' m 



~ 2.4 13G\. When 

'6 



and the one-loop QCD correction to the semileptonic decay gives /( 
we calculate the branching ratios of other rare processes, the branching ratio of B —>■ X^j 
which is experimentally measured with relatively high accuracy, must be considered as a 
prior constraint. 

In the SM, the CP asymmetry of i? ^ X^7 process 

T{B ^ X,7) - TiB ^ X 7) 



A^,(i?-^X7) 



is calculated to be rather small: A^.^ ~ 
js?! of the CP asymmetry, we have 



r(i?^x,7) + r(5^X7) ^^^^ 

0.5% 3^. By the recent experimental measurement 



0.30 < A^p{B X 7) < 0.14 



(22) 



at 95% C.L. In other word, the studies of the direct CP asymmetry in i? X^7 may uncover 
new sources of CP violation which lie outside the SM. Applying Eq. 1201 the CP asymmetry 
can be written as 



^^,(5-^X7) 



\CM\' 



^ - f{viz) + b{z,S)){l+^^)]lrn[CMCM)] 

ts tb 

l + ^^W,(/xJC;(^jlj(,23) 



v*v 



where z = {m^/mj^, v{z) and b{z,6) can be found in [3^ 



B. The decay width for two rare B processes 



1. l+l- {I = ijL, t) 



The decay constant of the pseudoscalar meson B^ is defined as [1 



{Q\-s^,^,h\BXP))=ip,f,^ . 



(24) 



With the equation of motion for quark fields, Eq. (|24|). one can write 

i2 



{Q\-s^,h\BXP))=if, 
13 



m,+m 

b ' s 



(25) 



Correspondingly, we derive the branching ratio as 

^2 q2 \y*y |2 ^ 4m^ 

Bs 

+ \F^+2m,FA, (26) 



with rrig and denote mass and hfe time of the meson respectively, and 



^ 2 BsJb,]^ m,+m 



^ 2 s.^s. I m, + m 

b ' s 

= —fsX^io - C[,) . (27) 
A point should be noted that there a CP asymmetry is observable in this process 

r(5 ^ r/) - r(fi ^ /I) 



2Im(e^jX 



(28) 



(i + ieopP)(i + ^f) ' 

and as indicated in the formula, it is induced by the mixing of 5^ and B^ |39[. Here, 

and Arrig^ is the mass difference in 5^ — B^ mixing, F^^ denotes decay width of the meson 
B^. 



2. B ^ Kl+l- {I = /i, r) 



According to Ref. |j4(|, the nonzero hadronic matrix elements for the exclusive decay 
B — > Kl^l^ are written as 



2 2 

ml — m~ 



{K{k)\s^,b\B{p)) = U{q'){2p-q),+ {h{q')-f^{q')) " , " 



(ir(A:)|w,.g^6|5(p)) = - ^^y^ {q\2p - q), - (mj - ml)q, 

2 2 
_ 777 — ?77 

{K{k)\-sb\B{p)) = i^foiq') , (30) 

m, — m 

b s 
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where = (p — k)'^ is the four-momentum transferred to the dilepton system. The resulting 
form factors are parameterized as 



/o(s) = /o(0) exp 



B B B 



his) = UiO) exp let-. + + c+(-^)= 



/t(^) = /T(0)exp 



SB B 



(31) 



With the effective Lagrangian Eq. (jl2j) and the hadronic matrix elements Eq. (j3U|) . we can 
write the transition matrix elements for B Kl^l^ as following 



(32) 



where p^ denotes the four-momentum of the initial B meson, and the form factors Fi are 



Lorentz- invariant. Following Ref. 



28| |. we define 9 as the angle between the three- momenta 



and in the center-of-mass frame of the dilepton. The energy-angular distribution of 
the decay products is 

dr(B^Kl + l-) a\Gl\V;Vj\l^ 2 2 N^rr^2,r.,2 ,r. 

^ - ^ " " ->^Hml,ml,s)P,{s[pf\F^\^+\F^ 



ds d cos 9 



2%5m| 



\F r + IF 



+4m^ml\Fj^ + 2m, {ml,m'^^, s)(3^Re{FgF*) cos9 
+ iml - ml + s)ReiF^F2)]} , 



(33) 



with s = {p^^ +Pi-f ^ A = V ^' c) = + fe^ + - 2{ab + be + ca). The 

kinematic quantities s, cos 9 have natural bounds 



4m^ < s < {nig — rrij^) , —1 < cos^ < 1 



(34) 



A particularly interesting quantity is the forward-backward asymmetry 

j' dcos9-r^^ - j" dcos9- '^^ 

n lis: d rns H — 1 C 



^^b(s) 



ds d cos ( 

It 



ds d cos 9 



f dcos9 , f . + f dco?,9-^ 

J ds d cos •> -1 ds d 



dV 



^2 p2 |T/*T/ |2 jp 



where the dilepton invariant mass spectrum ^ is given as 



(35) 



ds 



B 



-\^{ml,ml,s)(5\s[(5^^\Ff + \F, 
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2m^ 
1 + '- 



+ -X{ml,ml,s) 



+Am^m^ IFA + 2m, (m^ 

I b' A\ l\ B 



m 



I + s)ReiF,F*j} . 



(36) 



The form factors are formulated as 



-(m -m^)/o(s)^ 



b s 



m 



m 



Ms) -Ms] 



+ -{m - m)fo{s)^ ^ 



m, — m 

b s 



(37) 



rrig + rrij^ 

Our main interest is in the average forward-backward asymmetry < A^g >, which can be 
achieved from EqOni by integrating out the numerator and denominator separately over 
the dilepton invariant mass-square s. Basing on the preparations made above, we can 
numerically analyze the effects of various CP phases on the rare B decays in next section. 



IV. NUMERICAL RESULTS 

In this section, we present our numerical results for the rates of the rare B decays in 
the MSSM. As mentioned above, the most stringent constraint on the parameter space of 
'new physics' beyond the SM is the experimental bound on the branching ratio of i? ^ ^sl'- 

2. X 10"^ < BR{B X^7) < 4.5 x 10"^ 

In our numerical analysis, we take this as a constraint for the parameter space of the CP vio- 

I, 

lating MSSM. The inputs of the SM sector are j4^] a^^ = 1./128.8, = 80.23GeV, = 
91.18GeV, a^{m^) = 0.117, mf°''^ = 173.8GeV. The on-shell running masses of top and 
bottom quarks are related to the corresponding pole masses through j23| 



m,(mf''^) 



b\ (, / 



m 



pole 



5/(37r)a^(m- 



pole^ 



m 



pole 



5/(37r)a,(mP°'<=) ' 



(38) 
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FIG. 3: Taking tan/3 = 20, and /^^ = 300 (GeV), AM^_^ varies with c/i^p where (a)solid hne stands 
for (j)^p = 9^ as weh as 9.^ = 6^ = 0, (b)dash hne stands for (j)^p = 9^ as weh as 9^=9^= 0, and 
(c)dot hne stands for cp^p = 9^ as weh as 9^ = 9^ = 0, respectively. 
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FIG. 4: Taking tan/3 = 50, and /^^ = 100 (GeV), AM^^ varies with (\)^p where (a)solid hne stands 
for 4>^p = 9^ as weh as 9.^ = 9^ = 0, (b)dash hne stands for (j)^p = 9^ as weh as 9^=9^= 0, and 
(c)dot hne stands for cp^p = 9^ as weh as 9^ = 9^ = 0, respectively. 



For the physical CKM matrix elements, we adopt the Wolfenstein parametrization and the 
corresponding parameters are set as A = 0.85, A = 0.22, p = 0.21, 77 = 0.34 which stand 



4l[. In the hadronic sector, rrip^ 



for the central values permitted by present experiments 
5.37GeV, = 5.28GeV, = O.SOGeV, and r^^ = 1.46 x lO'^^^, = 1.54 x lO^^^^, 
and the decay constants are fp = fp^ = 0.21GeV. 

Beside the constraint from the branching ratio of BR{B -^s7), other strong constraints 
on the new CP violating phases originate from the — B^ mixing and resultant mass 
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FIG. 5: BR{B X^j) and A^p{B X^j) when mP"^" = 4.8GeV, tan/? = 20, as well as 
IJ,fj = 300 (GeV). In the figure, (a)solid line represents rigorous two loop analysis at (p^p = 6^ 
and 0^=0^= 0, (b)dash line represents one loop results plus threshold radiative corrections at 
^cp ~ ^2 ^3 = = 0) (c)dash-dot line represents rigorous two loop analysis at <p^p = 0^ 
and 62 = 0^ = 0, (d)dot line represents one loop results plus threshold radiative corrections at 
^cp ~ ^3 ^2 = = 0> (c)short-dash line represents rigorous two loop analysis at (p^p = 0^ and 
^2 = = 0, (f) dash-dot-dot line represents one loop results plus threshold radiative corrections at 
(pcp = 9, and 9^=6^= 0. 

difference: 

AM^^ > 9.48 X 10-^2 (GeV) 
as well as the mass difference in — mixing: 

AM^ = (3.304 ± 0.046) x 10"^^ (GeV). 
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With 1TL^^ = 300 (GeV) and tan/5 > 10, the theoretical estimation of the contributions to 
AMg from the SM and two Higgs-doublet sectors fits the experimental bound very well. 

d 

For the SUSY sector, the sfermion which contributes to AM^ , may belong to either the 
first generation or the third one. All the physical quantities which we are interested in, 
do not depend on the first generation parameters. Namely we are free to make our choice, 
no matter what parameters we set for the third generation of sfermions, we can assume 
that the new CP phases come from the first generation. Adding the SUSY contributions 
to that from SM and charged Higgs bosons, we still can make the total result satisfying 
the present experimental bound. Moreover, for the theoretical prediction on AM^ , we 
only consider the contributions of one-loop box diagrams and the threshold radiative cor- 
rections at large tan/3 2l|. Beside those parameters in the SM and two-Higgs-doublet 
sectors, the following supersymmetric parameters should be involved in our calculations 
u„, m„,, m. , 171- , m., rri-, rri-, as well as the trilinear Yukawa couplings A,, A,, A . 

A^H' 2,3' (32 3 -^2 3 * b ' ' f O t1 i, 5 s 

So far, there are no model-independent constraints on the masses of supersymmetric particle 
within the framework of CP conservative MSSM. Based on the following assumptions 

• Xi (7) is the lightest supersymmetric particle, 

• except t and 6, all scalar quarks are assumed to be degenerate in mass, i.e. = 

D q' 
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with the Tevatron data, loose bounds on the Stop and Sbottom masses are found 

m- > 140 GeV or < 64 GeV , 

m, > 210 GeV or < 32 GeV . (39) 

When we take into account the effect of the supersymmetric CP phases, those bounds would 
be further relaxed. In spite of this fact, we still take the bounds in Eq. ()39|1 seriously in our 
later numerical computations. Without losing too much generality, we fix the supersymmet- 
ric parameters as: \m^ \ = \m.j \ = 300 (GeV), = = 10 (TeV), m. = 1 (TeV), A^ = 
(GeV), m^^ = m.^ = 400 (GeV), m. = 500 (GeV), and m, = |4| = |AJ = 200 (GeV) 
unless otherwise noted. Now, our numerical results should be affected by the following CP 
violating phases: 9^ = arg(/i^), 6^ = arg(mj, 6^ = aTg{m^), 6^ = arg(v4J, 6^ = arg(v4J. 
In order to reduce degrees of freedom further, we set 9^=9^^ = 0. One of the reasons why 
we take this assumption on the parameter space is that the loop calculation for diagrams 
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inducing the lepton and neutron's EDMs restricts the argument to be < vr/ (5 tan j3) when 
those scalar fermions of the first generation are heavy enough. Additionally, we find that our 
numerical results depend on the argument 6^ rather moderately. Moreover, for the model we 
employ here, the mass of the lightest Higgs boson sets a strorig constraint on the parameter 
space of the new physics. As indicated in the literature lia, |2l|, the CP viola- 

tion would cause changes to the neutral-Higgs-quark coupling, neutral Higgs-gauge-boson 
coupling and self-coupling of Higgs boson. The present experimental lower bound for the 
mass of the lightest Higgs boson is relaxed to 60 GeV. In our numerical analysis we take this 
constraint for the parameter space into account. Now, we present our numerical results item 
by item. Since the present experimental result of AM^ constrains the 'new' CP phases in 
our calculation strongly, we discuss the mass difference in — mixing firstly. 

Taking tan/3 = 20, /i^ = 300 (GeV), we plot AM^ versus the CP phases (p^p in FiglHl 
where the solid line stands for (f)^p = 9^ and 0^ = 9^= 0, the dash line stands for (p^p = 6^ 
as well as ^2 — = O5 ^'^d ^^e dot line stands for = 9^ and 9^ =9.^ =0, respectively. 
FiglHis similar to FiglSl except there tan/5 = 50, /x^ = 100 (GeV). Under our assumptions 
on the supersymetric parameter space, the theoretical prediction on AM^ respects the 
experimental bound. 

For inclusive decay B ^^7; the present experimental observation on the branching 
ratio sets a constraint for the parameter space. Furthermore, the CP asymmetry in this 
process is highly sensitive to new CP violating phases because the SM contribution is only 
~ 0.5%. Taking tan/? = 20, /i^ = 300 (GeV), we plot BR{B X 7) and A^p{B X 7) 
versus the CP phases (p^p in FiglHl where the solid and dash lines represent (p^^p = 9^ and 

= 6'j = 0, the dash-dot and dot lines represent (f)^p = 9^ and 9^ = 9^ = 0, and the short- 
dash and dash-dot-dot lines represent (p^p = 9^ and 9^ =6*3 =0, respectively. From this 
figure, it is easy to note that there are very obvious differences between the exact two loop 
analysis (solid line for (p^p = 9^ , dash-dot line for (p^p = 9^ , and short-dash line for (p^p = 9^ , 
respectively) and the theoretical results which include one-loop contributions and threshold 
radiative corrections (dash line for (p^p = 9^ , dot line for (p^p = 9.^ , and dash-dot-dot line for 
(p^p = 9^, respectively). Taking (p^p = 9^ as an example, the CP asymmetry from rigorous 
two loop analysis can reach 6%, and that from threshold radiative corrections is smaller than 
1%, i.e. they are rather apart from each other, whereas the theoretical predictions made 
in the two scenarios on the branching ratios do not conflict with the present experimental 
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FIG. 6: BR{B X^j) and A^p{B X^j) when mP"^" = 4.8GeV, tanP = 50, as well as 
IJ,^ = 100 (GeV). In this figure, (a)solid line represents rigorous two loop analysis at = 9^ 
and 0^ = Ot = 0, (b)dash line represents one loop results plus threshold radiative corrections at 
^cp ~ ^2 ^i^d Oz — — 0, (c)dash-dot line represents rigorous two loop analysis at = 9^ 
and 9^=9^= 0, (d)dot line represents one loop results plus threshold radiative corrections at 
'Pcp ~ ^3 9.^ = 9^ = 0, (e)short-dash line represents rigorous two loop analysis at (j)^p = 9^ and 
^2 = 6*3 = 0, (f) dash-dot-dot line represents one loop results plus threshold radiative corrections at 
(t)cp = 9^ and 9^=9^ = 0. 

bound. 

CP violation is induced by both Standard Model (SM) sector and SUSY sector. In the SM, 
the CP violating parameter 77 = 0.34 (in the Wolfenstein parametrization, and corresponds 
to the central value permitted by the present data) which indeed induces non-zero CP 
violation in concerned processes. Therefore, even the SUSY phase takes special values as 
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FIG. 7: BR{B^ and A^-.p{B^ when mP°^'' = 4.8GeV, tan/3 = 20 as well as 

= 300 (GeV). In the figure (a)solid line stands for rigorous two loop analysis at (p^p = 9^ 
and 0^ = Ot = 0, (b)dasli line stands for one loop results plus threshold radiative corrections at 
^cp ~ ^2 ^J^d ^3 = 9^ =0, (c) dash-dot line stands for rigorous two loop analysis at (f)^p = 9^ and 
9^ = 9^ = 0, (d)dot line stands for one loop results plus threshold radiative corrections at = 9^ 
and 9 2 = 9^ = 0, (e) short-dash line stands stands for rigorous two loop analysis at 4>^,p = 6^ and 
9^ = 9.^ = 0, (f)dash-dot-dot line stands for one loop results plus threshold radiative corrections at 
= 9, and 9^=9^= 0. 

4>C'P = 0, TT, 27r, CP asymmetry may still exist. However, for the process B Xsj, the CP 
violation induced by the SM 1] is much smaller than 1%, by contraries, for the rare decay 
Bg —>■ the 77— induced CP asymmetry is greater than 1%, and the effect is expected 

to be observable even as the SUSY phase 0cp = 0, tt, 27r and does not contribute. In a 
word, the process B — > Xg^ provides a window for detecting the new CP sources beside that 
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existing in the CKM matrix. 

Similar to FigEl except for tan/3 = 50, = 100 (GeV), we plot BR{B X^'y) and 
A^p{B X^'-f) versus the CP phases (p^p in FiglHl With the setting for the parameter 
space, the branching ratio BR{B X^'-j) varies drastically with the CP phases = ^2*) 
but depends on the CP phase (p^p = 6^ very gently. At (p^p = 6^ = 7r/4, the CP asymmetry 
of rigorous two loop analysis can reach 5%, and that of threshold radiative corrections is 
about 3%, meanwhile the theoretical predictions on the branching ratios in both scenarios 
satisfy the present experimental bound. When cp^p = 0,^^, the CP asymmetry from the strict 
two loop analysis can reach 4%, and that from threshold radiative corrections is less than 
1%. 

Now, we present our numerical results on the rare decays B^ {I = fi, t). The 

present experimental upper bound on the branching ratio is BR{B^ f^^f^ ) < 1-5 x 10~^ 
at 90% C.L. Q- Taking tan/3 = 20, /i^ = 300 (GeV), we plot BR{B^ fi+fi-) and 
A^p{B^ /U^A* ) versus the CP phases (p^p in Fig|7l where the solid and dash lines stand 
for (p^p = 6^ and 9^ = 9^ = 0, the dash-dot and dot lines stand for (p^p = 9^ and 9^ = 9^ = 0, 
and the short-dash and dash-dot-dot lines stand for (p^p = 9^ and 9^ =9^ =0, respectively. 
From this figure, it is easy to find that there are evident differences between exact two loop 
analysis (solid line for (p^p = 9^ , dash-dot line for (p^p = 9^ , and short-dash line for (p^p = 9^ , 
respectively) and the theoretical results which originate from one-loop calculations plus 
threshold radiative corrections (dash line for (p^p = 9^, dot line for (p^p = 9^, and dash-dot- 
dot line for (p^p = 9^, respectively). For (p^p = 9^ ~ 37r/2, the theoretical prediction on the 
branching ratio BR{B^ — > fi^fi~) from one loop results plus threshold radiative corrections 
(dash line) can reach 1.2 x 10^® approximately, the exact two loop analysis (solid line) 
modifies the branching ratio to 5 x 10^^, and the CP asymmetry A^p{B^ /^^A*") is about 
3%. However, it is very difficult to detect this CP asymmetry A^p{B^ A*^/^") iii near 
future experiments with such a small branching ratio. As indicated above, the SM CP-odd 
parameter rj = 0.34 (in Wolfstein parametrization) induces a non-vanishing CP asymmetry 
when supersymmetric CP phases (p^p =0, tt, 27r. 

For larger tan/3, the situation is drastically different. Taking tan/3 = 50, /i^ = 
100 (GeV), we plot BR{B^ — > fi'^fi~) and A^p{B^ A^^A* ) versus the CP phases (p^p 
in FiglHl Clearly, there are obvious differences between the results of one loop contribution 
plus threshold radiative corrections and that of corresponding rigorous two loop calculations. 
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FIG. 8: BR{B^ and ^,,^(5^ ^ when m^"'^ = 4.8GeV, tan/3 = 50 as well as 

IJ,^ = 100 (GeV). In the figure, (a)solid line stands for rigorous two loop analysis at (p^p = 6^ 
and 0^=0^= 0, (b)dash line stands for one loop results plus threshold radiative corrections at 
^cp ~ ^2 ^3 = 6^ =0, (c) dash-dot line stands for rigorous two loop analysis at = 9^ and 
$2 = 0^ = 0, (d)dot line stands for one loop results plus threshold radiative corrections at (f)^p = 0^ 
and ^2 = 6^ =0, (e) short-dash line stands stands for rigorous two loop analysis at (p^^p = 6^ and 
9^ = 6,^ = 0, (f)dash-dot-dot line stands for one loop results plus threshold radiative corrections at 
(i)cp = 0, and 6*2 = 6*3 = 0. 

Additionally, the rigorous two loop prediction on the branching ratio BR{B^ /^^/W") sur- 
passes 10~^. Assuming the CP asymmetry to be induced by the CP phase (f)^p = 6^, the two 
loop result for the CP asymmetry A^p{B^ —>■ is about 5%, and the corresponding 

branching ratio is about 4 x 10~^. Although it is very challenging, this CP asymmetry could 
be hopefully measured in forthcoming experiments. 
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FIG. 9: Taking tan/3 = 20, /x^ = 300 (GeV), the branching ratio BR{B^ and CP 

asymmetry A^p{B^ vary with the CP phase (p^p = 9^, where (a)sond hne represents 

rigorous two- loop analysis with mP°^'^ = 4.6 (GeV), (b)dash line represents one-loop result plus 
threshold radiative correction with m^"^^ = 4.6 (GeV), (c)dash-dot line represents rigorous two- 
loop analysis with m^°'^ = 4.9 (GeV), (d)dot line represents one-loop result plus threshold radiative 
correction with mP°^^ = 4.9 (GeV). 

In Fig|7|and FiglHl we choose the pole mass of b-quark as m^"'*^ = 4.8 (GeV). Since the 
hadronic matrix elements depend on b-quark mass, we let vary within a certain range 
which is allowed by the data and see how it affects the branching ratio BR{B^ f^^t^~) 
and the CP asymmetry A^p{B^ fi^fi^). Taking tan/3 = 20, /i^ = 300 (GeV), we plot 
the branching ratio BR{B^ f^^f^~) and CP asymmetry A^p{B^ — * t^^f^ ) versus = 6^ 
in Figini To investigate the impaction of the b-quark mass on the measurable quantities, 
in the figure, we set it as m^"''^ = 4.6, and 4.9 (GeV), which correspond to the minimal 
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FIG. 10: BR{B^ r+r") and A^p{B^ r+r") when mP"^" = 4.8GeV, tan/3 = 20 as well as 
IJ,fj = 300 (GeV). In the figure, (a)solid line represents rigorous two loop analysis at (p^p = 6^ 
and 0^=0^= 0, (b)dash line represents one loop results plus threshold radiative corrections at 
^cp ~ ^2 ^3 = = 0) (c)dash-dot line represents rigorous two loop analysis at <p^p = 0^ 
and 62 = 0^ = 0, (d)dot line represents one loop results plus threshold radiative corrections at 
^cp ~ ^3 ^2 = = 0> (c)short-dash line represents rigorous two loop analysis at (p^p = 0^ and 
^2 = = 0, (f) dash-dot-dot line represents one loop results plus threshold radiative corrections at 
(pcp = 9, and 9^=6^= 0. 

and maximal values permitted by the present experiments respectively. Particularly, the 
branching ratio BR[B^ IJ'^ is relatively more sensitive to the b-quark mass. The 
present experimental errors approximately result in 5% theoretical uncertainty. As for the 
CP asymmetry A^p{B^ the corresponding theoretical uncertainty is only about 

2%. 
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FIG. 11: BR{B^ T+T-) and A^p{B^ r+r") when = 4.8GeV, tan/3 = 50 as well as 

fXjj = 100 (GeV) In the figure, (a)solid line represents rigorous two loop analysis at (j)^p = 9^ 
and 9^=9^= 0, (b)dash line represents one loop results plus threshold radiative corrections at 
(j)^p = 02 and 9 2 = 9^ =0, (c) dash-dot line represents rigorous two loop analysis at (p^p = 9^ 
and 02 — — 0) (d)dot line represents one loop results plus threshold radiative corrections at 
4'cp ~ ^3 ^^'^ ^2 — = 0) (e)short-dash line represents rigorous two loop analysis at (/)^p = 6^ and 
9 2 =6*3 =0, (f) dash-dot-dot line represents one loop results plus threshold radiative corrections at 
4'cp ~ ^i^d 9^ = 9^ = 0. 



Comparing with BR{B^ /i+yU^), the branching ratio BR{B^ t^t^) is enhanced 
strongly because r mass is much heavier than mass. Taking tan/5 = 20, = 300 (GeV), 
we plot BR{B^ r+r") and A^p{B_, T~^r~) versus the CP phases (p^p in FigHUl We find 
that the exact two loop predictions on the branching ratio BR{B^ r+r^) can reach ~ 10~^ 
approximately, and the CP asymmetry A^p{B^ t^t^) is about ±3% correspondingly. 
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When tan/? = 50, the branching ratio BR{B^ t^t^) is enhanced further. We plot 
BR{B^ — s> r+r^) and A^p{B^ t^t^) versus the CP phases (p^p = 6,^^^ in FiglTH with 
tan/5 = 50, fi^ = 100 (GeV). Assuming that the CP asymmetry is induced by the complex 
trilinear coupling A^, the two loop theoretical prediction on the branching ratio BR{B^ — >■ 
r+r~) is about 8 x 10^^, whereas A^p{B^ t^t^) ~ 5% at 9^ = 37r/2. Certainly, it 
is difficult to experimentally measure the rare decay B^ t^t^ . Similarly, the present 
experimental error for b-quark mass causes a theoretical uncertainties of ~ 5% for BR{B^ — >■ 
r+r^), and ~ 2% for A^p{B^ t^t^) respectively. 

We now discuss the branching ratio and forward-backward asymmetry in decays B — 
Kl^l~ {I = fi, t). The form factors for B —>■ Kl~^l~ (1 = /i, r) decays are given in Tablejl] 
which correspond to the central values presented in |40| . 



/o(0) 


/+(0) 


/t(0) 


^0 




T 


^2 




T 


^0 

C3 


C3 


T 


0.319 


0.319 


0.355 


0.633 


1.465 


1.478 


-0.095 


0.372 


0.373 


0.591 


0.782 


0.700 



TABLE I: The parameters for the parametrization Eql3H they correspond to the central values 
presented in Ref . |4C ] . 



Taking tan/5 = 50, fi„ = 100 (GeV), we plot BR{B ^ Kfi+fi-) and < Ap^ > {B ^ 
Kfi^fi~) versus the various CP phases (p^^p in Figll2| where the solid and dash lines stand for 
0CP ~ ^2 9^ = 9^ = 0, the dash-dot and dot lines stand for (p^p = 9.^ and 9.^ = 9^ = 0, and 
the short-dash and dash-dot-dot lines stand for cp^p = 9^ and 9^ = 9,^ = 0, respectively. This 
figure explicitly indicates that there are also obvious differences between the strict two loop 
results (solid line for (p^p = 9^, dash-dot line for cp^p = 9^, and short-dash line for (p^p = 9^, 
respectively) and the theoretical predictions which come from one-loop contributions plus 
threshold radiative corrections (dash line for (p^p = 9^, dot line for (p^p = 9,^, and dash- 
dot-dot line for (p^p = 9^, respectively). It is noted that although the branching ratio 
BR{B Kfi^fi~) can reach 3 x 10"'', the average forward-backward asymmetry is too 
small to be detected in any forthcoming experiment. 

Probably, the most interesting object to study is the rare decay B — » Kt^t~ because 
the average forward-backward asymmetry < App > {B ^ Kt^t^) is much larger than 
< App > {B ^ Kfi+fi-). Taking tan/5 = 20, fi^ = 300 (GeV), we plot BR{B Kt+t-) 
and < App > {B ^ Kt^t^) versus the CP phases (p^p in FigHSl Certainly, there are 
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FIG. 12: BR{B Kfi+iJ.-) and A^p{B Kfi+fi-) when = 4.8GeV, tan/3 = 50 as well 

as fi^ = 100 (GeV). In the figure, (a)sohd line stands for rigorous two loop analysis at (p^p = 6^ 
and 0^ = Ot = 0, (b)dash line stands for one loop results plus threshold radiative corrections at 
^cp ~ ^2 ^3 = 6^ =0, (c) dash-dot line stands for rigorous two loop analysis at (f)^p = 0^ and 
^2 = ^4 = 0, (d)dot line stands for one loop results plus threshold radiative corrections at (t)(~,p = 9^ 
and 02 = 0^ =0, (e) short-dash line stands stands for rigorous two loop analysis at 4>f^p = 9^ and 
9^ = 6.^ = 0, (f)dash-dot-dot line stands for one loop results plus threshold radiative corrections at 
= 9, and 9^=6^= 0. 

evident differences between the theoretical prediction of the exact two loop calculations 
and that of one loop result plus threshold radiative corrections. When (p^p = 9^ = 7r/2, 
if considering the one loop contributions plus threshold radiative corrections, one can have 
the branching ratio as BR{B — > Kt~^t~) ~ 10~^, whereas the strict two-loop calculations 
modify it to 8 x 10~^. The average forward-backward asymmetries < A^^ > (S — > Kt'^t~) 
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FIG. 13: BR{B Kt+t-) and A^p{B Kt+t') when mP°'^ = 4.8GeV, tan/? = 20 as well as 
= 300 (GeV). In this figure, (a)solid line stands for rigorous two loop analysis at (f)^p = 9^ 
and 9^=9^= 0, (b)dash line stands for one loop results plus threshold radiative corrections at 
(j)^p = 9^ and 9^ = 9^ = 0, (c)dash-dot line stands for rigorous two loop analysis at (j)^p = 9^ and 
^2 = 6t = Oi (d)dot line stands for one loop results plus threshold radiative corrections at cp^^p = 9^ 
and ^2 = = 0, (e) short-dash line stands stands for rigorous two loop analysis at (p^p = 9^ and 
^2 = ^3 =0, (f) dash-dot-dot line stands for one loop results plus threshold radiative corrections at 
^cp ~ 9^ = 9^ = 0. 

are about 10% for both cases. For tan f3 = 50, the branching ratio is enhanced further. When 
tan/5 = 50, fipi- = 100 (GeV) (Fig lT^ . the exact two loop calculations predict the branching 
ratio as BR{B —>■ Kt^t^) ~ 5 x 10~^, meanwhile the average forward-backward asymmetry 
< App > (B ^ Kt^t^) ~ 20%. If we can accumulate 10^° B mesons in experiments, one 
would be able to detect this asymmetry. 
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FIG. 14: BR{B Kt+t-) and A^p{B Kt+t') when mP°'^ = 4.8GeV, tan/? = 50 as well as 
= 100 (GeV). In this figure, (a)solid line stands for rigorous two loop analysis at (f)^p = 9^ 
and 9^=9^= 0, (b)dash line stands for one loop results plus threshold radiative corrections at 
(j)^p = 9^ and 9^ = 9^ = 0, (c)dash-dot line stands for rigorous two loop analysis at (j)^p = 9^ and 
^2 = 6t = Oi (d)dot line stands for one loop results plus threshold radiative corrections at cp^^p = 9^ 
and ^2 = = 0, (e) short-dash line stands stands for rigorous two loop analysis at (p^p = 9^ and 
^2 = ^3 =0, (f) dash-dot-dot line stands for one loop results plus threshold radiative corrections at 
^cp ~ 9^ = 9^ = 0. 

In Fig El and FiglT^ we take the pole mass of b-quark as m^"'*^ = 4.8 (GeV). The 
hadronic matrix elements depend on concrete values of b-quark mass. Taking tan j3 = 
20, fi^ = 300 (GeV), we plot the branching ratio BR{B Kt^t~) and the average 
forward-backward asymmetry < Apg > {B —>■ Kt^t^) versus (p^p = 6^ in Fig|T31 As for 
the b-quark mass, we set it as mP°^^ = 4.6 and 4.9 (GeV) respectively. Particularly, the 
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FIG. 15: Taking tan/3 = 20 as well as = 300 (GeV), the branching ratio BR{B Kt+t') 
and average forward-backward asymmetry A^gi^B Kt^t^) vary with the CP phase (j)^p = 6^, 
where (a)sohd hne represents rigorous two-loop analysis with m^"''^ = 4.6 (GeV); (b)dash line 
represents one-loop result plus threshold radiative correction with m^"'^ = 4.6 (GeV); (c)dash-dot 
line represents rigorous two- loop analysis with m^°'^ = 4.9 (GeV); (d)dot line represents one-loop 
result plus threshold radiative correction with m^"'^ = 4.9 (GeV). 

experimental uncertainty for the b-quark mass leads to a theoretical uncertainty of ~5% 
for the branching ratio BR{B — > Xr+r"), and ~2% for the average forward-backward 
asymmetry < A^^ > {B ^ Kt'^t'). 

The present experimental upper bound of BR{B^ — f^'^f^~) < 1.5 x 10"'' sets a strin- 
gent restriction on the supersymmetric parameter space. Along with the improvement of 
experimental precision and the accumulation of experimental data, this upper bound will 
be further modified. If so, we can expect that the new upper bound may lead to a concrete 
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FIG. 16: Taking mP^°''^ = 4.8GeV, 6*2 = ^3 = = Oj the correlation between /x^ and tan/3. 
In this figure, (a)the gray regions are permitted by the experimental bound on the branching 
ratio BR{B — > X^7), (b)solid line represents BR{B^ fJ-^fJ- ) = 5 x 10~®, dash line represents 
BR{B^ fi'^fi') = 10~^ and dot line represents BR{B^ M^/^") = 1-5 x 10"'^. 



constraint on the parameter space of our model. Using the exact two loop results, we plot 
the correlation of fi^ and tan/9 with the bound of the branching ratio BR{B ^sl) 
^2 = ^3 = = in Fig^l (a). The gray region is permitted by the present experiments. 
Corresponding to the two loop results for the branching ratios of BR{B^ — > we plot 

the correlation between /i^ and tan/? in Fib. El(b). Similar to Fig. IT^ except for 6^ = 7r/2, 
we plot the correlation of fi^ and tan/3 with the bound of the branching ratio BR{B X^j) 
in FiglTTI fa), and possible new upper bound on the branching ratio BR{B^ f^^f^ ) (Fig- 
ITTT b)). When 6^ = 0, the rare decay B^ fJ'^f^^ "will lead to a concrete constraint on the pa- 
rameter space if the new experimental upper bound on the branching ratio BR{B^ f^~^t^~) 
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FIG. 17: Taking m^"'^ = 4.8GeV, 9^ = 9., = 0, 9, = 7r/2, the correlation between fj,^ and 
tan p. In the figure, (a)the gray regions are permitted by the experimental bound on the branching 
ratio BR{B X^^), (b)solid line represents BR{B^ fi^fi^) = 5 x 10^®, dash line represents 
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reaches 10~^. As for the case 6^ = 7t/2, the branching ratio BR{B^ /^^A* ) 10~'' level 
will raise a stronger constraint on the model discussed here. Certainly, we plot Fig. and 
Fig. El under the hypothesis that the scalar quarks of the third generation are relatively 
light. If we push the scalar quark masses of the third generation to > ITeV, the situation 
would change drastically. Here we do not discuss such cases any further. 

Now let us simply discuss the dependance of the branching ratios and CP asymmetries 
(or forward-back asymmetries) on squark masses in the rare processes!?^ l^l^ and B — >■ 
Kl'^l' . Taking tan/5 = 40 and /i^ = —50 (GeV), we plot the branching ratio BR{B^ — »• 
fi^fj'^) as well as the CP asymmetry A^p{B^ f^^fJ' ) versus the mass of right handed 
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FIG. 18: Dependence of the branching ratio BR{B^ — > fi^fi^) and CP asymmetry A^p{B^ 
/i^/x^) on the right handed scalar quark mass m., with tan/3 = 40 and /x^ = —50 (GeV). In this 
figure, (a)solid hne stands for rigorous two loop analysis at 6^ = tt/2 and 9^ = 6^ = 0, (b)dash 
line stands for one loop results plus threshold radiative corrections at 62 = tt/2 and 9^ = 9^ = 0, 
(c)dash-dot line stands for rigorous two loop analysis at ^3 = 7r/2 and 9^=9^= 0, (d)dot line 
stands for one loop results plus threshold radiative corrections at 9^ = tt/I and 9,^=9^= 0, 

(e) short-dash line stands stands for rigorous two loop analysis at 9^ = Tr/2 and 9^=9^= 0, 

(f) dash-dot-dot line stands for one loop results plus threshold radiative corrections at 9^ = — 7r/2 
and 9^ = 9^ = 0. 

scalar top m. in Fig. ^1 Owing to the interference between the contributions of left 
handed and right handed stop, there is a resonant peak at m. = 400 GeV. When m. > 
1.5 TeV, the dependance of the branching ratio BR{B^ —>■ fi^fi^) and the CP asymmetry 
A^p{B^ — s> fJ'^fi^) on m. is very gentle. With the same choice of the parameter space. 
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we plot the branching ratio BR{B —>■ Kr^r^) as well as the Forward-Back asymmetry 
Apg{B Kt^t^) versus the mass of right handed scalar top m. in Fig. ^\ The resonance 
at m. = 400 GeV also originates from the interference between the contributions of left 
handed and right handed stop. For the strict two-loop theoretical predictions on BR{B — >■ 
Kt^t"), there is a small resonance around m. = 1 TeV which is due to the interference 
between the contributions of right handed stop and that of right handed scalar s-quark. A 
similar discussion about the dependence of the branching ratio BR{B X^'y) and the CP 
asymmetry A^p{B X^'j) on squark masses was given in our previous work [25[, and here 
we omit repetitions. It is noted that at large tan f3 scenarios, the corrections from sbottom 
on the rare process b sl^l^ only originate from two-loop box diagrams. This is the 
reason why the dependence of those branching ratios, CP asymmetries, and Forward-Back 
asymmetries on the sbottom masses is very weak. 



V. CONCLUSION 



Considering the constraints from the branching ratio BR{B — > X^j) and AM^^, we 
discuss the rare processes B^ — >■ /x^/i^, B^ — > r+r^, B Kfi^fi^ and B Kt^t~ in the 
CP violating MSSM at large tan (3. We find that there are evident differences between the 
theoretical predictions of exact two loop calculations and that of one loop contributions plus 
threshold radiative corrections. Additionally, the branching ratio BR{B^ —>■ fj^~^fi~) of exact 
two loop calculations can exceed 10~^, while the CP asymmetry can reach 3%. Although 
the two loop analysis predicts BR{B Kfi^fi^) > 10^'', the average forward-backward 
asymmetry < A^g > {B Kfi^fi^) is too small to be detected in future experiments. 
For the rare decay B^ t^t^, one has BR{B^ r+r^) > 10~^, meanwhile A^p{B^ — >■ 
r+r^) ~ 3%. Maybe, the most interesting object to study is the branching ratio BR{B — *■ 
Kt^t^) > 10^'^ along with a large average forward-backward asymmetry < A^^ > (i? — >■ 
Kt^t^) ~ 20%. As a by-product, we also find that the CP asymmetry of inclusive decay 
A^p[B —>■ X^j) can reach 5%, which is much larger than the SM prediction also. 
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FIG. 19: Dependence of the branching ratio BR{B — > Kt'^t") and average forward-backward 
asymmetry Apg{B Kt'^t'^) on the right handed scalar quark mass m-, with tan/? = 40 and 
= —50 (GeV). In this figure, (a)sohd hne stands for rigorous two loop analysis at 0^ = 7r/2 
and O3 = Ot — 0) (b)dash line stands for one loop results plus threshold radiative corrections at 
^2 = 7r/2 and ^3 = = 0, (c) dash-dot line stands for rigorous two loop analysis at ^3 = 7r/2 and 
02 = 0^ = 0, (d)dot line stands for one loop results plus threshold radiative corrections at 0^ = 7r/2 
and 6^=9^ = 0, (e)short-dash line stands stands for rigorous two loop analysis at 6^ = 7r/2 and 
9^ = 9^ = 0, (f)dash-dot-dot line stands for one loop results plus threshold radiative corrections at 
= -7r/2 and 9^=9^ = 0. 
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APPENDIX A: THE ONE-LOOP WILSON COEFFICIENTS AT THE WEAK 
SCALE 



The Wilson coefficients at the weak scale in Eq. ()12|1 are given as 
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The one-loop integrands are defined as 
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APPENDIX B: THE TWO-LOOP CORRECTIONS FROM THE Z, H PENGUIN- 
DIAGRAMS 
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APPENDIX C: THE TWO-LOOP CORRECTIONS FROM THE BOX DIA- 
GRAMS 
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Here, some two-loop functions are defined as 
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or the functions "if^^ as well as $(xo, Xi, X2) can be found in our forthcoming work 
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